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Abstract 



In strong (quasi-)Abelian fields, even at the one-loop level of the coupling 
constant, quantum fluctuations of fermions induce an effective Lagrangian 
density whose imaginary (absorptive) part is purely nonperturbative and 
known to be responsible for the fermion-antifermion pair creation. On the 
other hand, the induced real (dispersive) part has perturbative and nonper- 
turbative contributions. In the one-loop case, we argue how to separate the 
two contributions from each other for any strength of the field. We show nu- 
merically that the nonperturbative contributions are in general comparable 
with or larger than the induced perturbative ones. We arrive at qualita- 
tively similar conclusions also for the induced energy density. Further, we 
investigate numerically the quasianalytic continuation of the perturbative re- 
sults into the nonperturbative sector, by employing (modified) Borel-Pade. 
It turns out that in the case at hand, we have to integrate over renormalon 
singularities, but there is no renormalon ambiguity involved. 
PACS number(s): ll.15.Bt, lO.lO.Jj, ll.15.Tk, ll.80.Fv, 12.20.-m 

I. INTRODUCTION 

It has been well known for some time that the effects of the fermionic quantum fluctua- 
tions in space-time uniform Abelian gauge fields can be effectively integrated out, resulting 
in a one-loop effective action [|l|- 0. The results have been formulated also for the co- 
variantly homogeneous, thus quasi- Abelian, fields of the SU{2) gauge group and for 
specific nonhomogeneous magnetic field configurations 0. All the results can be extended 
to the case of the quantum fluctuations of scalar particles. The problems arising when gen- 
uinely non- Abelian fields with translationally invariant gauge-invariants are present [e.g., 
in S'f/(3)c] were discussed, e.g., in Refs. [p^pT|. 
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The quantum fluctuations of the strong gauge fleld itself (photons, or gluons) modify 
additionally those Lagrangian densities induced by the fermionic quantum fluctuations. Such 



two-loop effects have been successfully derived by Ritus [0 for homogeneous Abelian flelds, 
and further discussed by others [0,0] • In QED, such two-loop effects in the coupling 
constant change the one-loop result by at most a few per cent [ira). We will omit them in 
our investigation. 

There are basically two classes of phenomena associated with the presence of intense 
gauge flelds. 

Firstly, they can produce pairs of particles. For the Abelian case of strong homogeneous 
electric fleld, Sauter ||T5| showed this by investigating solutions of the Dirac equation in the 
corresponding potential,Q and Schwinger [§] by using methods of action integral. Green's 
functions and proper time. Differential probabilities for pair creation were investigated in 
Refs. [0 and |T8|. In the latter Reference, the quasi-Abelian model was applied to inves- 



tigation of the quark pair production in chromoelectric flux tubes. Experimental evidence 



related with the pair production in a strong QED (laser) fleld was reported in Ref. [0 . 

The pair production has its origin in the imaginary (absorptive) part of the effective 
Lagrangian induced by the fermionic quantum fluctuations in the strong fleld. That part is 
entirely nonperturbative in nature, because the production effects are ~exp(— const. /t?) and 
thus cannot be expanded in positive powers of the fleld-to-fermion coupling parameter g. 

On the other hand, the other class of phenomena is associated with the real (dispersive) 
part of the induced effective Lagrangian. In QED, this class includes the following phe- 
nomena that affect a low energy {u -C rUe) photon wave entering the region of the strong 
background fleld: photon splitting, change of the photon speed, and birefringence. Works 



on the theoretical aspects of these phenomena include Refs. 0]- |24]. The experimental 



aspects of birefringence in strong magnetic flelds are discussed in [^- |^. The dispersive 
part of the induced action leads in principle to those corrections of the classical Maxwell 
equations which originate from the (fermionic) quantum fluctuations. 

The aim of the present paper, while dealing with the dispersive part of the induced 
action, is somewhat different from these works. We concentrate on the concept of sepa- 
rating the nonperturbative from the perturbative contributions in the induced dispersive 
action when the product of the (quasi)electric fleld S and the coupling constant g is large: 
gS/rn^ ~ 1, where m is the fermion mass. Subsequently, we numerically investigate the 
two contributions. Afterwards, we use the discussed quantities as a "laboratory" for testing 
and investigating the efficiency of methods of quasianalytic continuation. The latter meth- 
ods, involving the (modified) Borel-Pade approximants, allow us to obtain approximately 
the nonperturbative contributions from the approximate knowledge of the perturbative con- 
tributions and by employing the Cauchy principal value prescription in the inverse Borel 
transformation (Laplace-Borel integral). These considerations can give us insights into 
the problems of extraction of nonperturbative physics from the knowledge of perturbative 
physics in gauge theories, in particular in various versions of QCD (high-fiavor, low-fiavor). 

In Section II, we argue how to perform the mentioned separation into the perturbative 



^ A related problem was first considered even earUer by Klein [0 who investigated solutions of 
the Dirac equation with a high vertical barrier potential (Klein's paradox). 
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and nonperturbative contributions. After identifying the two contributions, we investigate 
numerically their values for various values of the field parameter d ~ g£/m?. In Section 
III we then carry out an analogous analysis for the induced energy density, the latter being 
in principle observable. In Section IV we then numerically investigate, for the induced 
Lagrangian and energy densities, (quasi) analytic continuation from the perturbative into the 
nonperturbative sectors, employing the method of Borel-Pade for the induced Lagrangian 
and a modified Borel-Pade for the induced energy density. We encounter integration over 
renormalon poles, whose origin is nonperturbative, and we show how to carry it out. Section 
V summarizes our results and conclusions. 



II. INDUCED DISPERSIVE LAGRANGIAN DENSITY 



We start by considering the Euler-Heisenberg Lagrangian density [0 which is the 
part induced by the quantum fluctuations of the fermions in an arbitrarily strong (quasi- 
)Abelian homogeneous field 



5C 



— exp[— is(m^ — ie' 
s 



g^abcoth.{ags) cot{bgs) — ^{a^ — b^) — 



Here, g is the field-to-fermion coupling parameter (in electromagnetism it is the positron 
charge Cq), m is the mass of the (lightest) fermion, and the parameters a and b are Lorentz- 
invariant expressions characterizing the (quasi)electric and the (quasi)magnetic fields £ and 
B, respectively 



-8^ + ^ 



1/2 



1/2 



/V2. 
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We note that ab= \S-B\, and a'^—b^ = S^ — B^. Further, a— >■ when \B\ —>-0, and b—* \B\ when 
|£| — >-0. In the Lorentz frame where B\\S, we simply have: a= and 6= Expression (|I]) 
can be obtained, for example, directly by integrating out the fermionic degreees of freedom 
in the path integral expression of the full effective action, and employing the proper-time 
integral representation for the difference of logarithms. 

As denoted, the integration in (^ is performed along the positive real axis infinitesimally 
below it, avoiding in this way the poles on the real axis appearing due to the (Lorentz- 
invariant version of the) magnetic field b (function cot{bgs) there). If the path in (0) were 
above the real axis, then we would obtain a nonzero imaginary part of the Lagrangian density 
even in the pure magnetic field case. This would imply particle creation in this case, which 
is physically unacceptable. The path in (|l|) reproduces for the case of the pure magnetic 
field the known real density, and for the case of the pure electric field the known complex 
density 0. The path in (|I|) is suggested also from the extension of the formal approach of 
Ref. to the general ab^O case. Namely, when ab^O, we need to evaluate in this approach 
two traces: one trace originating from a^O and discussed in [their Eq. (4-116)]; and the 



other trace of the evolution operator of an harmonic oscillator, originating from 67^0 of 
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the form X]exp[— 215(76(71 + 1/2)] (n = 0, 1, . . .; s > 0). The latter trace becomes convergent 
after regularization s i— > s — ie (e— s>+0), i.e., the path in (|l|). 

Performing a contour integration in the fourth quadrant of the complex proper-time 
s-plane (cf. Fig. |ID, expression (H) can be rewritten as 



6C 



dz 



[z+ie] 



^2 

g^abcot {ag{z+ie)) coth {bg{z + ie)) H (a^- 

3 



{z + iey 



(4) 



where s = —iz + e now runs along the negative imaginary axis (and As seen from 

Fig. H the result (|D is actually independent of the precise path of the proper time variable 
s in (|lD, as long as s runs from near the origin towards s = +oo and passes each positive 
pole on the right, i.e. precisely the class of paths satisfying the physical condition that the 
pure magnetic fields cannot produce particles. In (|l|) and (|), the familiar counterterm 
oc {a^—b"^) [= {S'^—B'^)] is included which makes the integral finite. This divergent term leads 
to the renormalization of the field in the leading Lagrangian density C^'^^ = {£'^—B'^)/2. We 
now divide the integration region into intervals for the integration variable agz: iq = [0, vr/2], 
ii = [71/2, 37r/2], . . ., in = [(n — l/2)7r, (n+l/2)7r], . . . Each interval, except io, contains in its 
middle one pole of the integrand. The corresponding series for the real (dispersive) part of 
the Lagrangian density is 



ReSC = ReSCo + ^ ReSC 

n=l 

— exp 
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pcot{w) coth.{p{w+n7r)) 



pcot{w) coth(pw) H — (1— p 

a / I 

dw exp — — 
-7r/2 V a J [w+m:) 

p coth(pn7r) ^ (1— 



1 



(5) 
(6) 
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w 



f-7r/2 



+Re / dwexp — — 

J-7T/2 V a J [w + ie) 



3{w+mT) (w+nir)^ 

pcoth.{pmT) 1 



(n>l) . 



(7) 



Here, e' = ega +0, and we used the notation 



p= - - - 



a 



a 



6C = 6C/ 



' 4-2 ^ 

ma 



87r2 



(8) 



and we introduced the dimensionless integration variable w = agz when agz is in the interval 
io, and w = agz—n7i when agz is in the interval in (n> 1). In (|^), we separated the integrand 
into a part that is entirely nonsingular in the integration region, and a part that is singular 
but gives a finite value of integration since the Cauchy principal {V) value has to be taken. 
From a formal point of view, we note that 6Cq doesn't "feel" the poles of the integrand as 
depicted in Fig. |l|, while 5£„ {n>l) "feels" the pole s = —imT/{ag) via the principal value 
part in (0) that is proportional to 
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r'^l^ , exp(— w/a) 
dw— — — 

-7r/2 (W + ie'j 



= 7^ / — exp(— ti?/a) 
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-El 



TT 

25 



-Ei 



Ta 



-2 [x + xV(3! 3)+a;V(5! 5) + - ■ ■] 1^=^(25) if 5 > 1 
-(eVa;) + + - ■ ■] U=v(2a) if a < 1 
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We note that the dispersive part of the induced Lagrangian density as normahzed here 
(i)^® depends only on two dimensionless parameters - on parameter p = h/a which char- 
acterizes in a Lorentz-invariant manner the ratio of the strengths of the (quasi)electric and 
(quasi) magnetic fields [cf. and on parameter d = {ga)/m? which characterizes the 

combined strengths of the (quasi)electric field parameter a and the field-to-fermion coupling 
g. In the perturbative weak-field limit, a is small. In this case, when reintroducing in 
z = w/{ag) 



Re5£o 



— expl — 2;m 

z 



g^ 

g^ahcoi{agz) coi}i{hgz) + — (a^ — 6^) 



(10) 



we can see that the real part of expression (^) is approximately reproduced, since formally 
7r/{2ag) — oo. In this case, the conventional perturbative expansion of the dispersive La- 
grangian density in powers of g^ (i.e., inverse powers of x) can be performed (cf. [0], 0) 



^£pert. = (^cil\d^ + CsS! d^ + CgS! d^ + 

where the expansion coefficients are 



(11) 



C3 



945 



C5 



14175 



3{l-py + 22p\l-py + 19p 



etc. 



(12) 



Expression (|n|) can be derived alternatively by purely perturbative methods - the terms 
~a^" can be obtained by calculating the one-fermion-loop Feynman diagram with 2n photon 
external legs of zero momenta. Expression ([TT|) is a divergent asymptotic series and it gives 
the usual perturbative corrections to the Maxwell equations. On the other hand, the formal 
small-a expansion of Re6Co of (H) [or equivalently: (0)] reproduces the terms (|ll|) and yields 
in addition the terms ~aexp(— const./a). The latter terms may in principle be dangerous 
for the interpretation of Re5£o of @ as the perturbative part of the induced density, since 
they are nonanalytic and could thus signal physical nonperturbative effects. However, in 
the Appendix we demonstrate that these terms are only an artifact of the abruptness of 



the infrared (IR) proper-time cutoff z< 1/Afj^ (^ir = (2/7r)m a 



.0 These terms are 



^ The energy cutoff AiR~m\/S is low in the case when the perturbative effects dominate (i.e., at 
a < 1), but is higher when the nonperturbative effects are significant (at a > 1). The nonperturbative 
effects here reside in the infrared (IR) sector of (fermionic) momenta q < Air, and the effective 
contributing size of this sector gets larger when d grows. 
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therefore not of a physical nonperturbative origin. In the Appendix we further show that 
Re5£o of (iD should be reinterpreted as the limit with an infinitesimally softened IR cutoff, 
the latter limit being numerically the same but having no nonanalytic terms in the small-a 
expansion. That expansion is then identical to (pH]) . 

On the other hand, the densities ReSCn {n>l) of (|^ represent the nonperturbative part 
of the induced dispersive density @, for two reasons: 

• The integration over the proper-time z runs here in the vicinity (in fact, across) the 
n'th pole of the integrand of (p. The poles of the integrand are in the nonperturbative 
regions. We recall that these poles are also the source of the nonzero imaginary 
(absorptive) part of the density leading to the fermion-antifermion pair creation, a 
clearly nonperturbative phenomenon. 

• The densities Re5£„ {n>l), independently of the pole structure of their integrands, 
become appreciable only in the strong-field (large-a) regime while in the weak-field 
(small-a) regime they decrease faster than any power of a, i.e., they do not contribute 
to the perturbative series (pi]). Each of the two integrals in the curly brackets of 
(|^ behaves as ~ aexp[7r/(2a)] when a — > +0, and thus the entire Iie6Cn behaves as 
~aexp[— (n — l/2)7r/a] {n>l) in this limit.0 

Therefore, the nonperturbative effects contained in Re(5£„ (^>1) are of two types, one 
type being characterized by the poles of the integrand, and the other type by what we may 
call strong-field effects. From the above discussion, it further follows that we have some 
freedom in choosing the IR proper-time cutoff: 2;<1/Ajj^ is such that Ajj^ ~m^a and that 
all the possible poles must lie at z^s above the cutoff 1/Ajj^. We took Ajj^ = nm?a with 
fi; = 2/7r, but any n satisfying l/7r<K~l would be acceptable as well. 

From a somewhat different perspective, we can imagine transforming a truncated per- 
turbation expansion for Re(5£P''''*' /a of (|TTp (with several terms) via the Borel-Pade approx- 
imation. The resulting integrand approximately reproduces the integrand of (^), including 
the poles structure. Thus the integration over the ra'th pole, contained in Re5£„ of (|^), can 
be interpreted as the n'th renormalon in the density, i.e., a nonperturbative quantity. We 
will return to this point later in this paper. 

Thus, the densities (H) and (|^ result in the fermion-induced perturbative and non- 
perturbative contributions, respectively, to the Maxwell equations. The fields were taken, 
strictly speaking, to be homogeneous in space and time. In practical terms, this means that 
they are not allowed to change significantly on the distance and time scales of the Compton 
wavelength of the fermion 1/m. For electro-magnetic fields, m is the electron mass, and 
1/m is about 4 ■ 10~^^ m, and 1.3 • 10~^^ s. 

An indication of the relative size of the perturbative and nonperturbative fermion- 
induced corrections to the Maxwell equations can be obtained by comparing the correspond- 
ing contributions to the induced Lagrangian density. This is done in Figs. Figures || (a), 
(b) show the dimensionless perturbative (||) and nonperturbative (|^) induced Lagrangian 



^ If we did not take the principal Cauchy value in (^, but some other prescription (which in the 
case at hand would be wrong), Re5£„ would behave as ~exp[— n7r/a]. 
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densities, respectively, as functions of the (quasi)electric field parameter d (||), at four dif- 
ferent fixed values p = b/d of the magnetic-to-electric field ratio. The case of the pure 
(quasi) magnetic field (p.m.f.) is also included in the Figures, as function of b. For the 
p.m.f. case, we normalized the Lagrangian density in analogy with (^, i.e., S£ is obtained 
in that case by dividing 6C by m^6^/(87r^). The separation between the perturbative and 
the nonperturbative part was performed in the p.m.f. case analogously, i.e., the proper- 
time z <TT / {2bg) contributions were defined to be perturbative, and those from z>7T/{2bg) 
nonperturbative. We point out, however, that in the latter case the nonperturbative con- 
tributions do not involve the renormalon-type ("pole-type") effects, but only strong-field 
effects (cf. previous discussion). In Fig. the corresponding ratios of the nonperturbative 
and perturbative induced densities are presented.0 When moving beyond the perturbative 
region (i.e., when 1), we see from these Figures that the nonperturbative parts in general 
become relatively significant and often even dominant. 

Once we come into the nonperturbative regime (a~l), however, we must keep in mind 
that the pair creation, originating from the large absorptive part, will become so strong 
as to render the solutions of the corrected Maxwell equations unstable. We will quantify 
somewhat this fact in the next Section in the case of the induced energy density in QED. 

In Figs. (a), (b), the densities were normalized according to (|^), so that the tree-level 
reference values for the densities are 

^'°'-':'°v(^)=^(i-P^). (13) 

Therefore, increasing only the coupling parameter g, while leaving the (quasi) electric field a 
unchanged, results in correspondingly larger relative corrections originating from the induced 
parts, both nonperturbative and perturbative. In the special case of QED, on the other hand, 
g = eo is small [a = e^/ (An) 1/137], and the overall induced Lagrangian density accounts 
usually for less than 0.5 per mille of the total Lagrangian density when d<l (see the next 
Section on related points). 



III. INDUCED ENERGY DENSITY 

In this Section, we discuss the induced energy densities. Energy density is in principle a 
measurable quantity. It is not Lorentz-invariant. If the (quasi) electric and (quasi)magnetic 
fields are mutually parallel, the various induced energy densities can be obtained directly 
from the corresponding induced Lagrangian densities 



SUk = a 
6lAk = d 



da 
9Re(5£fc 



da 



+ Re5Ck (A; = 0,1,2,...) 



(14) 



^ The total induced dispersive Lagrangian densities, and values of the truncated perturbation 
series ( |Tl| ) (including ~a®), are included in Figs. ^ in Section IV. 



where we denoted, in analogy with (|) 



' A~1 ^ 

ma 



(15) 



With the restriction to parallel fields £ \\ B (i.e, \£\ = a and \B\ = b) we do not lose the 
generality since, for any configuration of S and B, there always exists a Lorentz boost, 
perpendicular to the plane of the fields, so that in the boosted frame the two fields are 
parallel. The corresponding perturbative and nonperturbative parts of the energy densities 
in such frames are 



Re5Wn 



^/"^ dw 



w 

— exp I — — 

w \ a 



^-1 
a 



1 1 
pcot(w) cothfpw) + -(1—p^) 



+ 



pcot{w) coth(pw) + p — 



w cot(w) ^ 2 
sinh^(pw) 3 



nrr 
exp I — — 



7T/2 / ^ 

ctwexp — — 
\ a 



(w+mr) 



a 



(16) 

pcot{w) coth.{p{w+mi)) 



(w+nTr) 



p coth(pn7r) (1— p^) 



w niT 

f-7r/2 



3{w+n'K) {w+mr)'^ 
[ dwexp (—-r) 

J-n/2 \ a J 

p ( coth(pn7r) 



, , / cothfnfty + riTr)) 
pcot(w) + 



P 



+ 



p 



nn 



+ 



p 



sinh (pmr) 

r coth(pr;,7r) H 

a smh [pmT) 

H — — coth(pn7r) sinh f — ^ I , 
niT \2aJ ) 



2 

+ 



W + nTTj 

1 



sinh {p{w + mi)) ^ 
2 



3 (w+mr) (w+mr)^ 













(25). 







(n > 1] 



The tree-level density in the normalization convention used [cf. (|^)] is 

47r2 



r 4-2 ^ 

ma 



9' 



The perturbative power expansion of the induced energy density 6U is 

^^pert. = (d_^i\a^^ d0 ~a^ + d^h\d^ + - 
where the expansion coefficients are 



, 1 

di = — 

45 

d^ = 



3 + 5j9^ - 



945 



10 + 21p^ - + 2/ 



1 



14175 



21 + 50p^ - 2V + 10/ - 3p^ 



etc. 



(17) 



(18) 



(19) 



(20) 



The results for the induced perturbative ([T6| ) and nonperturbative parts (p]?]), and their 
ratios, are presented in Figs. § (a)-(b) and ^, respectively, in analogy with Figs. |] (a)-(b) 
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and 1^. The case of the pure (quasi)magnetic field is not included in Figs. because in 
this case 5U = —6C and thus the relevant information is already contained in Figs. |^-^. The 
behavior of the induced energy densities is, in broad qualitative terms, similar to that of the 
induced Lagrangian densities.^ 

In the special case of QED, similarly as for the Lagrangian densities in the previous 
Section, the total induced energy densities account for a very small part of the total energy 
density (0.2-0.3 per mille when d^l) and can become significant only when the field becomes 
exceedingly large (a~10^). The same is true also for the dielectric permeability tensor Sij: 
In the direction of the fields, we have 6e\\ =e\\ — l = ad(Re6C)/da\h, i.e., by (|T^)-(|T5|) we 
have = {6U + R.e6C)a/{2iT), which is about 10"'^ for a ^ 1 and p = 1. Therefore, the 
effective coupling parameter along the field direction ay =a/e\\ changes by about one per 
mille, while a± = a/e± remains unchanged since e± = 1. Therefore, in QED, any quantity 
which can be expanded in powers of the coupling parameter alone (without fields) remains 
a perturbative quantity. QED then remains a perturbative theory despite such strong fields 



- cf. also Ref. ||28| on that point. 

The energy density is not stable in time when d ^ 0, due to the energy losses to pair 
creation of fermions of mass m. It decreases by about 50 percent in the time ti/2 



ti/2 ~ — exp +- 



p-K coth(p7r) 



- , (21) 

m 



where a = g'^/{ATT). The factor in the square brackets, appearing due to the presence of the 
(quasi) magnetic field, is usually not essential in the estimates since it is ~1 for p<5. In the 
case of QED and with p = 0, ti/2 is about 0.9 ■ 10^m~^, 0.4 ■ 10^m~^ and 0.3 ■ W^m~^ for 
a = 0.5, 1, and 5, respectively. Here, ^ 1.3 ■ 10"^^ s is the electron Compton time. 



IV. QUASIANALYTIC CONTINUATION INTO THE NONPERTURBATIVE 

SECTOR 

In this Section, we use the discussed induced densities as an example on which to test and 
get some insights into methods of approximate analytic (i.e., quasianalytic) continuation. In 
various physical contexts, such methods allow one to extract all or part of the information 
on the nonperturbative sector from the knowledge of the perturbative sector alone. We will 
use the method of Borel-Pade transformation, or a modification thereof. 

One may ask whether the perturbation expansions (|TT]) and ( [T^ ) allow us to obtain the 
full, including the nonperturbative, information about the corresponding densities. The 
answer for the Lagrangian density is yes, but under the condition that we take in the 
corresponding Borel-Pade approximants the Cauchy principal values when integrating over 
the positive poles of the Fade integrand in the inverse Borel transformation. This is refiected 
in the terms ie in the denominators of the integrands of (^ and/or (0). More specifically, 
we first Borel-transform (B) the perturbation series (pH]) 



^ The total induced energy densities, and values of the truncated perturbation series (19) that 
include terms ~a^, are included in Figs. ^ in Section IV. 



B 



Ci{p)a + C3{p)a +C5{p)a + 



(22) 



then construct an [N/M]B{a;p) Pade approximant to B of (|2^) ,P1 and then apply the inverse 
Borel transformation 



^ p[N/M] r^£pert 



{a;p) 



dw exp 



w 



[N/M]Biw;p) . 



(23) 



On the other hand, the real part of the actual density (^) can also be written as a Borel-type 
integral, when introducing w = agz and e' = age in and normalizing the density according 

to m 



Ke5L = Re / dw exp 
Jo 



w 
d 



w 



pcos[w 



smw + ie 



— coth(pw) + -(1-p 



(24) 



The expansion of the integrand of (|^ 
with the Borel transform (E2D with d 



excluding the exponential, in powers of w is identical 

^ ^ ^w, as it should be. Comparing ( p3D with the exact 

result (24), we see that the Borel-Pade method (|23|) will be efficient in (quasi) analytic con- 
tinuation if Pade approximants [N/M]b{w) approach the integrand of (0) in an increasingly 
wide integration interval of w when the Pade order indices and M {^N) increase. This 
in fact happens, since the integrand in ( P^ ) is a meromorphic function in the complex plane 
whose poles structure on the positive axis is especially simple - there are only single (not 
multiple) poles, located at w = 7i, 27r, Sn. Pade approximants to power expansions of such 
functions are known to approximate such functions increasingly better when the Pade order 
indices N^M increase |2^. Near the poles w^nn the integrand behaves as ~(w— riTr+ie')"^- 
Hence, for obtaining the real (dispersive) part of the density, the Borel integration over the 
poles has to be taken with the Cauchy principal value (CPV) prescription - not just in 
the exact expression (^), but also in the approximate expression (PB]). Thus, the Borel 
integration in (pSj) over the n'th pole, i.e., the n'th renormalon contribution, has in the case 
at hand no renormalon ambiguity. As the Pade order indices N^M are increased, we thus 
systematically approach the exact lieSC via the CPV of p3|). This means that in the case 



at hand [strong (quasi-) Abelian fields with fermionic fiuctuations included], the full induced 
Lagrangian density can be obtained on the basis of the knowledge of perturbation expansion 
(pA|) for weak fields and the CPV prescription. The more terms in ([TT| ) [and thus in (p2|)] 
we know, the higher Pade order indices N^M we can have, and hence the closer to the full 



Lagrangian density we can come via (pSj). 



On the basis of the knowledge of the first four nonzero perturbation terms in (|Tl]) and 
correspondingly in (p2|), we can construct the following Pade approximants of the pertur- 
bative Borel transform (|2^): [1/2]b, [1/4]b, [3/4]b- Then we can calculate the corresponding 



^ [N/M]{a) Pade approximant to (^) is defined by two properties: 1. it is a ratio of two polynomi- 
als in o, the nominator polynomial having the highest power d''^ and the denominator d^^; 2. when 
expanded in powers of d, it reproduces the coefficients at the terms d" in (22) for n<A^-|-M; it is 
based solely on the knowledge of these latter coefficients c„ {n<N+M). 
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Borel-Pade transforms via p3| ) with the CPV prescription. The corresponding results of the 
approximants for the full induced density R.e6C are presented in Fig. ^, together with the 
exact numerical values calculated by (p^ ) in Section II. The curves are given as functions of 
the (quasi)electric strength parameter d for four fixed values of the magnetic-to-electric field 
ratio p, and Fig. ^ (d) is for the case of the pure (quasi)magnetic field (a = 0). We see that 
the highest order ([3/4]) Pade-Borel results agree well with the exact results over the entire 
depicted region of a. When the Fade order indices N and M (~A^) increase, the region 
of agreement includes increasingly large values of d. For comparison, we also included the 
results of the truncated perturbation series (TPS) made up of the first four nonzero terms 
of ( pAD [in Fig. (d): for the corresponding p.m.f. case], i.e., those perturbation terms which 
the presented Borel-Pade transforms are based on. 

If we apply the very same procedure in the case of the energy density - Borel- 
transforming the series 6U^'^^^'/d of (|19]), constructing Fade approximants, and carrying 
out the inverse Borel transformation by using the Cauchy principal value (CPV) prescrip- 
tion - the results are disappointing. It turns out that increasing the Fade order indices N 
and M (~A^) does not generally result in a better precision. For example, for p^ 0.5 and 

a ^0.5, the Borel-Pade transforms of the order [3/4] and [3/6] give significantly worse results 
than those of the lower order [1/4]. The reason for this erratic behavior of the Borel-Pade 
approximants in this case lies in the more complicated poles structure of the Borel-Pade 
transforms. This can be seen if we rewrite 6U in the Borel-integral form analogous to (p4|), 
obtained from (|2^ ) by applying relation (|1^ 



ReShl = Re / dwexp 
Jo 



w 
d 



-I] 



w 



1 



1 



pw , , , J- ,^ n, 

'—^7 — coth(pw) + -{1+p ) + — 
sm [w+ie ) S 



(25) 



The expansion in powers of w of the integrand in (P3|), excluding the exponential, gives of 
course the exact Borel transform of the perturbation series (|1^) divided by d (and replacing 



1 

w 



pw , / N 1 /-, 2 

— — coth(pw) + -[l-p 
sm (w + ie) 3 



di{p)w + d3{p)w^ + d^{p)w^ + 



B 



1 

w"^ 

sw'^'^iw^p) 

w 



(26) 



However, we now see that this integrand has a double poles structure on the positive w axis, 
the double poles located at w = 7i, 27i, 37i, .... The Fade approximants to the power series 
(p6|) have great trouble simulating this double poles structure adequately. When they do 
it by creating one single or two nearby real poles, say near w = n, then it turns out that 
the inverse Borel transformation via the Cauchy principal value (CPV) prescription often 
gives good results. However, when the Fade approximants try to simulate the double pole 
near w = n hj creating two mutually complex-conjugate poles a± ib (a ^ tt, |6| <^ tt), the 
inverse Borel transformation gives very unsatisfactory results. This occurs, for example, in 
Fade approximants [3/4]{w;p) and [3/6]{w]p) for p<0.5. Heuristically we can understand 
that such a simulation is bad, because the structure of the integrand in ( pSj ) suggests that a 
double pole at a—ib alone, just below the real axis, would do a better job, but it is not allowed 
in the Fade approximants. The latter is true because the perturbation expansion (|26|) is 
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explicitly real for real w's, and this property is hence shared also by the Pade approximants, 
enforcing for each complex pole another pole which is complex-conjugate. 

To overcome this problem, the idea is to modify the Borel transformation of the pertur- 
bation series ([19|) in such a way that the resulting transformed series is represented by a 
(meromorphic) function without any double poles on the real positive axis, in contrast to 
the Borel-transformed series (^61) . This, in fact, can be implemented in the easiest way by 
using the following modification of the Borel transformation (MB): 



dMB 



^^pert. 



dw 

MB ' 



{w;p) 

= B 



w 



di{p)w + ds{p)w^ + d^{p)w^ + 



2 4 fi 

~ 1 111 in in 

5W^^'] {w;p) = d,{p)\ + + d,{p)— + ■■■ . (27) 



This trick changes every double pole in the B into the corresponding single pole in the MB. 
Then we apply Pade approximants [A^/M]mb('U^) to the MB series (0), and carry out the 
corresponding inverse modified Borel transformation 

MBP[^/^] \6W''''-] (fi; p) = - r dwe^^ (--\ \NIM]ub{w- p) (28) 

L J fi JO V fi / 

with the Cauchy principal value (CPV) prescription when integrating over the single poles. 
This CPV prescription originates again from the ie' terms in the double poles structure 
(if— nvr+ie')"^ of the Borel-transform (B) integrand of ( p5D that is now changed to the single 
poles structure (if — ravr + ie')"^ in the modified Borel-transform (MB) of (^Tf ) whose Pade 
approximants [N / M]mb{w;p) appear in (pSj). 

The numerics clearly confirm that these MBP's (^) are well behaved, i.e., they approx- 
imate well the actual full induced energy density SU{d;p) in the region of fi which is getting 
wider when the Pade order indices N and M {^N) increase. The results are presented in 
Fig. 1^, where the MBP's for the first three possible Pade order indices [2/2], [2/4] and [4/4], 
along with the exact numerical results, are shown as functions of fi, at four fixed values of 
p = b/d. Another reason why the results now behave better than those of the usual BP 
tranforms lies in the fact that the Pade approximants ([2/2], [2/4] and [4/4]) are now more 
diagonal than earlier ([1/2], [1/4] and [3/4]). This is due to one additional power of w in 
the MB series (P7|), as compared with the usual B series. The diagonal and near-diagonal 



Pade approximants are known to behave better than the (far) off-diagonal ones [^. In fact. 
Figs. 1^ suggest that clear improvement - extension of the fi range of agreement with the 
exact results - sets in when we switch from [2/2] to [4/4] MBP, while the off-diagonal [2/4] 
MBP may even be slightly worse than [2/2]. For comparison, we also included the results 
of the truncated perturbation series (TPS) made up of the first four nonzero terms (up to 
~fi®) of (|l9l), i.e., the terms on which the presented Borel-Pade transforms are based. The 
case of the pure (quasi) magnetic field was not included in these Figures because in this case 
6U = —6C and thus the information on this case is contained in Fig. |^ (d). 

This application of Borel-Pade transformations and their modification may give us some 
insights into how the (quasi) analytic continuation from the perturbative (small fi) into the 
nonperturbative (large fi) regions can be carried out in other theories whose exact behavior 
in the latter region is still theoretically unknown. One such example is the perturbative QCD 
(pQCD), where some observables are known at the next-to-next-to-leading order (N^LO). 
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The coupling parameter in that case [d i— > as{Q'^)] can be quite large when the relevant 
energies of the process are low (Q ~ 1 GeV), thus rendering the direct evaluation of the 
N^LO TPS unreasonable or at best unreliable. When applying Borel-Pade transformations 
or modifications thereof to such series, we are faced with two major problems: 

• The first problem is of a more technical nature. Since only very few, at most two, 
coefficients beyond the leading order are known, the Fade approximants associated 
with the (modified) Borel transform of the series have low order indices (A^, M < 2) 
and thus do not necessarily reproduce the location of the leading poles on the positive 
axis, if they exist, adequately. 

• The second problem is of a deeper theoretical nature. Knowing too little about the 
behavior of QCD in, or close to, the nonperturbative regime, we do not know how to 
integrate over the possible positive poles in the inverse (modified) Borel transformation 
- this can be termed the infrared renormalon ambiguity pOl . 



In the discussed case of integrated fermionic fiuctuations in strong (quasi-)Abelian fields 
- for the Lagrangian and energy densities - we do not face any of the two afore-mentioned 
problems since the exact solution is known. We have to apply the Cauchy principal value 
(CPV) prescription in the integration of the Borel-Pade transform of the induced dispersive 
Lagrangian density, and in the modified Borel-Pade transform of the induced energy density. 
The CPV is the direct consequence of the path (e parameter) in the exact solution (P [-^ @)]. 
The knowledge of the full theoretical solution in the latter case also tells us that the poles 
structure of the usual Borel transform of the induced energy density is more complicated 
(double poles), so that we have to apply a modified Borel-Pade transform which changes 
the double poles into a single poles structure. 

We point out that the positive poles - renormalons - discussed in the present work cannot 
be directly identified with the usual infrared (ultraviolet) renormalons in QCD (QED). The 
latter renormalons, as defined in the literature [^, are interpreted in the perturbative 
language as originating from renormalon chains at low (high) momenta k. The renormalon 
chains are momentum-/c gluon (photon) propagators with n chained one-loop insertions, 
where n can be arbitrarily large. In the model at hand, however, only quantum fiuctuations 
of fermions, in the slowly-varying strong fields, are considered; the effects of the quantum 
fluctuations of propagating gluons (photons) were not included in the discussed effective 
model. The positive poles, i.e. renormalons, in the present model originate from a collective 
effect of arbitrarily many very soft gluons (photons) coupling to a fermion loop or to a 



fermion propagator - cf. ||3T|. The relevant parameter of the effective coupling of these 
soft gauge bosons to the fermions, appearing in the induced effective action, is a = ga/rn? 
and it can be large due to the strong fleld a and/or due to the strong coupling g. These 
nonperturbative contributions are then roughly ~exp(— const. /a) = exp[— const. m^/((7a)] - 
cf. (p3|), (|28|). This is similar, but not identical, to the infrared renormalon contributions 
in QCD ~exp(— const. Y^f^). We may be tempted to term the renormalons discussed in the 
present paper as infrared renormalons due to their nonperturbative origin in the infrared, 
although this name is reserved for the afore-mentioned QCD-type renormalons. 

Various QCD and QED apphcations of the Borel-Pade approach, with CPV prescription, 
have been made in [0- The new method of Ref. ||3^ gives modifled real and imaginary 



parts of the Borel-Pade of 5C, in comparison to the usual CPV prescription, when the Pade 
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approximants [N/M]b have poles off tlie positive real axis. This may influence the speed 
of the convergence of the Borel-Pade transforms towards the full solution when the Pade 
order indices and M A^) increase. This method came to our attention after finishing 
the manuscript. 

Two other references [BSl- [BBl are also somewhat related to our work. Dunne and Hall 



35| considered, among other things, the question of resummation of the (one-loop) Euler- 



Heisenberg (EH) Lagrangian density by using the knowledge of the perturbation expansion 
of the Borel transform. Since they did not use Fade in addition, they needed at least an 
approximate information on all the coefficients of the series to reconstruct approximately the 
nonperturbative sector. Jentschura et al. on the other hand, did not employ the Borel 
transform, but applied directly to truncated perturbation series (TPS) of the EH density 
a numerical method (Weniger sequence transformation) which differs from Pade in several 
aspects. Their results of resummation are better than the direct application of Pade to the 
TPS of the EH Lagrangian, but they are worse than the results of the combined Borel-Pade 
method. 



V. CONCLUSIONS 

We introduced the concept of separation of the induced dispersive action into the non- 
perturbative and perturbative parts. We then investigated numerically the nonperturbative 
contributions to the dispersive (real) part of the Lagrangian density and to the real energy 
density, induced by quantum fluctuations of fermions in the strong (quasi-)Abelian fields 
that don't change significantly in space-time over the typical fermionic Compton wave- 
lengths 1/m. There are only nonperturbative contributions in the absorptive (imaginary) 
part of the strong field Lagrangian density, the latter part being responsible for the fermion- 
antifermion pair creation. On the other hand, the nonperturbative contributions in the real 
(dispersive) sector are in general also significant and can often even dominate over the per- 
turbative induced contributions there. The induced dispersive Lagrangian density modifies 
the Maxwell equations for strong fields. The induced energy density is in principle an ob- 
servable quantity. When the (quasi)electric fields are strong, however, these densities decay 
fast (in ~10'^ Compton times, for a~l). These two induced densities lead to a change in the 
dielectric permeabihty tensor of the vacuum. In the special case of QED, all these induced 
effects are below one per cent unless the fields are huge (a~10^). 

We then used the discussed induced quantities as a "laboratory" to test and investigate 
the efficiency of specific methods of quasianalytic continuation from the perturbative region 
(weak fields) into the nonperturbative region (strong fields). We employed the method of 
Borel-Pade for the induced dispersive Lagrangian density, since the function represented 
by the Borel transform series has only simple poles. For the induced energy density, we 
had to employ a modified Borel-Pade transformation since the function represented by the 
(nonmodified) Borel transform series has double poles. We found out numerically that 
such quasianalytic continuations become precise over an increasing region of the effective 
expansion parameter a when the number of available terms in the perturbative expansion 
increases. This means that the quasianalytic continuation gradually becomes the analytic 
(exact) continuation when the number of the perturbative expansion terms accounted for 
increases. The Borel integration over positive poles (renormalons) is necessary. The correct 
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prescription for the integration over these poles, in the case at hand, is the simplest one - 
the Cauchy principal value (CPV) prescription, its origin being the path (e parameter) in 
the exact solution (|I]) (^)]. Such analyses could give us some insight into the problems 
faced in QCD when nonperturbative contributions to observables are investigated either on 
the basis of the perturbative results themselves or by using other models that are at 
least partly motivated by perturbative methods. 

The correct analytic continuation, in the discussed case of strong background gauge 
fields, is the one employing the simplest (CPV) prescription for integration over the poles in 
the Laplace-Borel integral. This appears to be in agreement with the conclusions of Ref. 



which were obtained from quite different considerations involving the renormalization group 
- that the vacuum polarization induced by the intense gauge fields is in principle determined 
by the information on the behavior of the theory in the perturbative region. The situation 
in QCD is less clear. A necessary condition for the existence of the (correct) analytic con- 
tinuation from the perturbative into the nonperturbative regime in QCD is that a nontrivial 
infrared stable fixed point exist for the running strong coupling parameter. Such an infrared 
stable fixed point, however, seems to exist only if the number of the quark flavors is high 
{Nf > 9) ||3^. For the real (low-A^/) QCD, a phase transition takes place, and methods of 
analytic continuation have probably only a limited range of applicability. Stated differently, 
in this case the full knowledge of the perturbative sector probably does not allow us to 
obtain information on the deep nonperturbative sector. In such a case, it is probable that 
even the renormalon ambiguity in the low-flavor perturbative QCD (pQCD) is an intrinsic 
ambiguity that cannot be entirely eliminated with pQCD-related methods alone. 
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Appendix A. ON THE ANALYTICITY OF 6Co 

In this Appendix we will clarify the nature of the nonanalytic terms ~exp(— const./a) 
that appear in the naive expansion of the perturbative part Iie6CQ{d;p) of (P) around the 
point a = 0. Such terms may in principle be dangerous for our interpretation of (|^) as the 
perturbative part of the induced Lagrangian density, because they have the nonanalytic 
structure similar to those terms that appear in the nonperturbative parts Re(5£„(a;p) of 
(0), the latter containing genuinely nonperturbative contributions due to the singular (pole) 
structure of the integrand. We will show that the mentioned terms in @ are an artifact 
of having the abrupt infrared (IR) cutoff there, and that they disappear as soon as the 
abruptness of the infrared cutoff is (infinitesimally) softened. 

In the proper-time formalism, the IR and UV regions correspond to the high and the 
low values of the proper time, respectively [|^. In the proper-time integral (^) for SC, 
the IR region of large proper time z [z > n/{ag)] contains poles, the latter leading to 
nonperturbative effects. The region of smaller z has no such singularities and thus no 
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nonperturbative effects. Therefore, the perturbative part of 5C should cover the latter 
region, and suppress the IR region. The general way to do this is to introduce, in the spirit 
of approaches of a nonnegative regulator Pe(if) {w = agz) in the proper-time integral 



1- 



, (A.l) 



YleSCf = —Re / — Peiw) exp ( — — ) pcot(w + ie') coth(pw) H — (1— 
Jo w \ a J I 3 

where the minimal IR regularization requirements are 

Pe{w) ^ 1 for W <^ 1 , Pe{w) ^1 for W TT . (A. 2) 

The nonnegative parameter e indicates that we can choose a class of such regulators. In fact, 
we will require that for small e a large chunk of the perturbative region, namely the w-region 
of approximately [0,7r/2], survive in (|A.1|) . Thus we restrict the minimal conditions ( [A.2|) 



to the following ones, when e <^ 1: 

Pe{w)^l for w ^ 7r/2 - ^/e , Pe(w) < 1 for w ^ n/2 + ^/e . (A. 3) 

A seeming alternative to ( |A.1| ) would be to introduce a regulator Pei^z) that would scale 
as a function of z = w/{ag) instead of w. But this possibility must be discarded because 
then the condition of suppressing the pole structure [pe{z) <^ 1 for z > Tx/{ga) =7r/(m^a)] 
cannot be reconciled with the condition of the survival of a large chunk of the perturbative 
region [pe{,z) f» 1 for z < Ti/{2rn?a)\ at various values of a simultaneously. 

The conditions ( |A.3|) are designed in such a way that the limit e— *+0 would apparently 



lead to the abrupt IR regulator appearing in YieSCQia^p) of (^, with the abrupt cutoff at 
w = 7r/2. We can choose the following specific one-parameter family of regulators Peiw) 
satisfying the afore-mentioned conditions: 

, , Peiw) - / X 1 1 (w — -n/2\ 

Pe{w) = , Pe{w =--- arctan . A.4 

p,{0) 2 TT \ e J 

When e — * +0 (^t^O), these regulators differ from the abrupt cutoff regulator outside the 
narrow w-interval [tt/2 — -y/e, 7r/2 + ^/£] by at most r^^/€ 

r 1 - ieM{n/2 - w)-' + 0{e') if «; < f -^^ , 1 
P'^"^^ - \ {e/7r){w - 7r/2)-i + 0{e^) ifw>l + ^,j ^^'^^ 

while they may differ from the abrupt version significantly only in the afore-mentioned 
narrow interval. The first thing to check would be that the regularized expression ( |A.1|) 



with the regulator ( |A.4|) -( |^l5| ), in the limit (^7^0) really gives numerically the result 

(I) of the abrupt cutoff. Stated otherwise, we should check that the limg^+o in front of the 
integral ( |A.1| ) can be moved into the integral, without changing the result. For such a check, 
we need to see that the contributions in ( [A.l|) from the singular (poles) regions {w > vr/2) 



are suppressed toward zero when Such a check is straightforward and we performed 

it. It turns out that the w-regions [(n— l/2)7r, (?7,+l/2)7r] around the ra'th pole Wn = nn are 
suppressed by a factor when n > 2, and by at least a factor ^^/e when n = l. Thus, all 
these contributions go to zero when On the other hand, on the ly-interval [0,7r/2], 

there are no singularities of the integrand and the regulator is virtually equal to 1 in the 
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entire interval when e +0. Therefore, on this interval we can automatically push the 
limiting procedure into the integral. Thus we really have 



lim ReS£.''f'\a]p) = Re6Co{a]p) 



(A.6) 



i.e., the numerical value of the perturbative part with the infinitesimally "softened" IR cutoff 
is the same as that of the perturbative part with the abrupt IR cutoff (|^). 

Now we will investigate the expansions of the above two expressions around the point 
5 = 0, in order to see the difference in the (non)analiticity structure between the two cases. 
We can find the small-a expansion of Re5£o(S;p) of by expanding first the integrand 
(without the exponent) there, i.e., the Borel transform, in powers of w. As argued in Section 
IV [cf. Eqs. (p2D-(plD], this expansion yields ( p2D with d \—>- w, where Cj(p)'s are given by 
(|l^). Then the term-by-term integration over w leads to the small-a expansion of Re(5£o 



Re5£o(a;p)^'"''-^ 
= 'ci{p)l\d^ + cs{p)3\d^ 
+0 (a^exp[-7r/(2a)] 



r7T/2 ^ r-K/2 

ci{p) I dw e'x.'^{—w / a)w + C2,{p) I dwexp{—w/a)w + 
Jo Jo 

3 



aexp 

^ * 2d 



TT 



(0 +C3b) (^) + 



(A.7) 



Incidentally, the coefficient at aexp[— 7r/(2a)], written as an infinite sum, is just the value 
of the Borel transform at w = tt/2 [cf. remark following Eq. (pi))] 



Clip) ( +C3{p) 



+ 



2x 2 



(A.8) 



The coefficients of terms (9 (a^ exp[— 7r/(2a)]) can be obtained in an analogous manner, by 
using derivatives of the Borel transform with respect to w at w = 7r/2. Expressions ( A.7 )- 
( |A.8| ) show explicitly the following: In the small-a expansion of Re6Co of (|), in addition to 
the usual perturbation expansion part ([TT|) that is analytic at a = 0, we obtain formally also 
terms ~a" exp[— vr/ (2a)] which are nonanalytic at a = 0. One might suspect that such terms 
could possibly be of nonperturbative origin, and below we will show that they are not. More 
specifically, we will show that they are an artifact of the abruptness of the IR cutoff and 
that they are de facto not there, in the sense that they disappear when we consider instead 
of Re(5£o its numerical equivalent, i.e., the e —>■ +0 limit of the left-hand side of (|A.6|) . To 
show this, we have to expand the latter expression [at e^O - i.e. ( [A.!] )] around a = 0. For 
that, we first Taylor-expand the regulator pe{w) ([A.4|) , which is analytic everywhere|] in 
powers of w for small e 



n 



w - — 

2 Vtt 



n+2 



(A.9) 
(A.IO) 



'''in contrast to the abrupt cutoff when pq{w) = 1 for w<{tt/2), and po(w)=0 for w>{it/2). 
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The other part of the integrand in ([A.l|) , without the exponent, is the Borel transform whose 
small-a expansion is ( p2|) with a ^ w. Combining this and (|A.9|) -( [A.10 ), we obtain after 
some straightforward algebra^ the small-a expansion of ( |A.1| ) around a = for small e 



(exp.) 



;;6 



ci{p)l\a +C3(p)3!a +c^{p)h\a 



+ 
+ 



- Ci(p)+C3(p) 

vr 



TT 



ci(p) 



2x2 

vr 



4! + 

C3(P) +C5(p) 



5!a^ 



6! + ■ ■ ■ + 



(A.ll) 



Here we see explicitly that the small-a expansion of 'K'&buf'^ (a; p) of ( |A.1|) at nonzero e exists 
and that this function is analytic there, having no nonanalytic terms ~exp(— const./a), in 
contrast to the expansion of Re5£o where e was set equal to zero exactly (i.e., inside the 
integral). Further, expansion ( [A. lip goes over into the usual perturbation expansion bCP^'^'^' 
(|n|) when 

These considerations thus lead us to the following conclusions: 

• The perturbative part of the induced Lagrangian density, Re5£o as defined in (^, 
has an abrupt IR cutoff at w = 7r/2, and it is numerically equal to the corresponding 
expression with an infinitesimally softened IR cutoff - cf. left-hand side of ( A. 6 ). 



The small-a expansion of Re5£o(o;p) reproduces the usual perturbation expansion 
(|TT]) plus nonanalytic terms ~a" exp (—const. /a) [cf. ( [A.7| )]. 



The small-a expansion of the corresponding expression ( |A.1|) with a softened IR cutoff 
[e 7^ 0) yields no nonanalytic terms; when the softening of the IR cutoff becomes 
infinitesimal (e — +0, £7^0), the expansion becomes identical with that of the usual 
perturbation expansion (pT]). 



The above points show that the nonanalytic terms in the small-a expansion of 
Re5£o(o;p) are only an artifact of the abruptness of the IR cutoff (the cutoff reg- 
ulator becomes a nonanalytic function of the proper time vj) and are thus not of a 
nonperturbative physical origin. Re5£o(0';p) should be reinterpreted as the limit with 
the infinitesimally softened IR cutoff [the left-hand side of ( |A.6| )], the latter being 
numerically the same but its small-a expansion having no nonanalytic terms. 



We again integrate term-by-term; and we repeatedly use the identity: /g°° duexp{—u)u'^ = n\. 

18 



REFERENCES 



[1] W. Heisenberg and H. Euler, Z. Physik 98, 714 (1936); 

[2] V. Weisskopf, Kgl. Danske Videnskab. Selskahs. Mat.-fys. Medd. 14, No. 6 (1936). 
[3] J. Schwinger, Phys. Rev. 82, 664 (1951). 

[4] W. Greiner, B. Miiller, and J. Rafelski, Quantum electrodynamics of strong fields, 594 

pp. (Berlin, Heidelberg, Springer 1985); 
[5] W. Dittrich and M. Renter, Effective Lagrangians in Quantum Electrodynamics, 244 

pp. (Springer 1985). 

[6] J. Schwinger, Particles, Sources and Fields, Vol. II, 306 pp.. Chapters 4-8 (Addison- 
Wesley, 1989). 

[7] C. Itzykson and J.-B. Zuber, Quantum Field Theory, 705 pp.. Chapter 4-3 (McGraw- 
Hill, New York, 1980). 

[8] I. A. Batalin, S. G. Matinyan, and G. K. Savvidy, 5*0^;. J. Nucl. Phys. 26, No. 2, 214 
(1977). 

[9] R. Ragazzon, Phys. Rev. D52, 2422 (1995). 
[10] R Schwab, Phys. Lett. B109, 47 (1982). 

[11] A. B. Migdal and S. B. Khokhlachev, Pis'ma Zh. Eksp. Teor. Fiz. 41, No. 4, 159 (1985) 
[JETP Lett. 41, No. 4, 194 (1985)]. 

[12] V. I. Ritus, Zh. Eksp. Teor. Fzz. 69, 1517 (1975) [Sov. Phys. JETP 42, 774 (1976)]; V. I. 
Ritus, Zh. Eksp. Teor. Fiz. 73, 807 (1977) [Sov. Phys. JETP 46, 423 (1977)]; V. I. Ritus, 
Effective Lagrange function of intense electromagnetic field in QED, in Proceedings of 
the conference Frontier tests of QED and physics of the vacuum, pp. 11-28, Sandansky, 
Bulgaria, June 1998, eds. E. Zavattini, D. Bakalov and C. Rizzo (Heron Press, Sofia, 
1998), also |hep-th/9812T24 . 

[13] M. Renter, M. G. Schmidt, and C. Schubert, Ann. Phys. (N.Y.) 259, 313 (1997). 

[14] G. V. Dunne and C. Schubert, Nucl. Phys. B564, 591 (2000) 

[15] F. Sauter, Z. Phys. 69, 742 (1931); ihid 73, 547 (1931). 

[16] O. Klein, Z. Phys. 53, 157 (1929). 

[17] A. I. Nikishov, Zh. Eksp. Teor. Fiz. 57, 1210 (1969) [Sov. Phys. JETP 30, 660 (1970)]. 

[18] A. Casher, H. Neuberger, and S. Nussinov, Phys. Rev. D20, 179 (1979). 

[19] D. L. Burke et al, Phys. Rev. Lett. 79, 1626 (1997). The electron pair production 
measurement reported in this work included, beside several soft ( "homogeneous" ) laser 
photons, one energetic Compton-backscattered photon. The latter photon thus probes 
the field inhomogeneity since its wavelength is much smaller than that of the laser field. 

[20] Z. Bialynicka-Birula and I. Bialynicki-Birula, Phys. Rev. D2, 2341 (1970). 

[21] S. L. Adler, Ann. Phys. (N.Y.) 67, 599 (1971). 

[22] E. Brezin and C. Itzykson, Phys. Rev. D3, 618 (1971). 

[23] W.-Y. Tsai and T. Erber, Phys. Rev. DIO, 492 (1974); zMd 12, 1132 (1975). 

[24] W. Dittrich and H. Gies, Phys. Rev. D58, 025004 (1998); W. Dittrich and H. Gies, 
Vacuum birefringence in strong magnetic fields, in Proceedings of the conference Fron- 
tier tests of QED and physics of the vacuum, pp. 29-43, Sandansky, Bulgaria, June 



1998, eds. E. Zavattini, D. Bakalov and C. Rizzo (Heron Press, Sofia, 1998), also |hep- 
th/980641?l. 



[25] D. Bakalov et al, PVLAS: Vacuum birefringence and production and detection of nearly 



19 



massless, weakly coupled particles by optical techniques, Nucl. Phys. B Proc. Suppl. 35, 
180 (1994). 

[26] S. A. Lee et ai, Measurement of the magnetically-induced QED birefringence of 
the vacuum and an improved laboratory search for light pseudoscalars, FERMILAB- 
PROPOSAL-P-877A, 1998. 

[27] G. Cvetic, private notes. 

[28] R. D. Peccei, J. Sola, and C. Wetterich, Phys. Rev. D37, 2492 (1988). 

[29] George A. Baker, Jr. and Peter Graves-Morris, Pade Approximants , 2nd edition, 746 

pp. (Encyclopedia of Mathematics and Its Applications, Vol. 59), edited by Gian-Carlo 

Rota (Cambridge University Press, 1996). 
[30] A. H. Mueller, Nucl. Phys. B250, 327 (1985); V. I. Zakharov, Nucl. Phys. B385, 452 

(1992); A. I. Vainshtein and V. I. Zakharov, Phys. Rev. Lett. 73, 1207 (1994); zbzd 75, 

3588E (1995); M. Beneke, Phys. Rept. 317, 1 (1999) and references therein. 
[31] C. B. Chiu and S. Nussinov, Phys. Rev. D20, 945 (1979), cf. especially Fig. 1 there. 
[32] P. A. R§czka, Phys. Rev. D43, R9 (1991). 
[33] M. Pindor, |hep-ph/9903T5T . 

[34] U. D. Jentschura, Phys. Rev. D62, 076001 (2000). 

[35] G. V. Dunne and T. M. Hall, Phys. Rev. D60, 065002 (1999). 

[36] U. D. Jentschura, J. Becher, E. J. Weniger, and G. SofF, |hep-ph/9911265 . 

[37] M. Beneke and V. M. Braun, Phys. Lett. B348, 513 (1995); P. Ball, M. Beneke and 

V. M. Braun, Nucl. Phys. B452, 563 (1995); Y. L. Dokshitzer, G. Marchesini, and B. 

R. Webber, Nucl. Phys. B469, 93 (1996); E. Gardi and G. Grunberg, JHEP 11, 016 

(1999). 

[38] S. G. Matinyan and G. K. Savvidy, Nucl. Phys. B134, 539 (1978). 

[39] E. Gardi, G. Grunberg, and M. Karliner, JHEP 07, 007 (1998); F. A. Chishtie, V. Elias, 

V. A. Miransky, and T. G. Steele, |hep-ph/9905291 . 
[40] R. D. Ball, Phys. Rep. 182, 1 (1989); G. Cvetic, Ann. Phys. (N.Y.) 255, 165 (1997). 



20 



(s-plane) 



ag 



-2i 



ag 



-3i 



ag 



FIG. 1. The contour integration, in the complex s-plane, needed to rewrite in the form (^. The 
location of the poles is denoted explicitly. 
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FIG. 2. (a) Perturbative and (b) nonperturbative induced dispersive (Euler-Heisenberg) Lagrangian 
densities [cf. (^) and (0)] as functions of the (quasi)electric field parameter a (||), at various fixed values of 
the magnetic-to-electric field ratio p = b/d (^). The actual values of the curves for p«0, p= 1.0 and p = 5.0 
have been multiplied here by factors 10, 5 and 1/2, respectively, for better visibility. Included is also the 
case of the pure (quasi)magnetic field (p.m.f.), for which the x-axis represents b — gb/m^. 
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FIG. 3. Ratios of the nonperturbative and the perturbative induced dispersive Lagrangian densities for 
the cases depicted in Figs. 0. For the p.m.f. case, the x-axis represents b = gb/ni?' . 
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FIG. 4. (a) Perturbative and (b) nonperturbative induced energy densities [cf. ( p^ ) and (|T^] as func- 
tions of a at various fixed values of p = 6/a (||). The actual values of the curves have been multiplied, for 
better visibility, by the denoted factors, just as in Figs. ^. 




FIG. 5. Ratios of the nonperturbative and the perturbative induced energy densities for the cases 
depicted in Figs. ^. The ratio for p — 1 varies strongly for d — 0.5-1.5 because the perturbative induced 
density has a zero at a w 1 . 1 . 
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FIG. 6. Borel-Pade approximants (BP's) to the induced dispersive Lagrangian density (| 
of a, for various values of p = 6/a: (a) p = 0; (b) p = 0.5; (c) p = 5.0; (d) pure magnetic field (a = 



) as functions 
0). Depicted 

are those BP's (23) which are based on the Pade approximants [1/2], [1/4] and [3/4] of the Borel-transform 
(^^. The numerically exact curves [sum of curves of Figs. H] (a) and (b)] are also included and they virtually 
agree with the [3/4] BP results. Included are also the the results of the truncated perturbative series which 
include terms [in (d): ^b^]. 
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FIG. 7. Modified Borel-Pade approximants [MBP's - cf. (p8|)] to the induced energy densities, based 
on the Pade approximants [2/2], [2/4] and [4/4] for the MBP's (p7|), as functions of a, at fixed values of 
p = b/a: (a) p = 0; (b) p = 0.5; (c) p=1.0; (d) p = 5.0. 
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